load("pairData.Rda")

##############################################
## unidimensional pairwise comparison model fixed alpha
##############################################
library(MCMCpack)
set.seed(481113)

out.1d.fixedalpha <- MCMCpaircompare(pairData, 
                                     theta.constraints=list('1015' = 0.25,
                                                            '1014' = "-"),
                                     alpha.fixed=TRUE,
                                     a=0, A=1,
                                     verbose=1000,
                                     burnin=20000, mcmc=100000, thin=10,
                                     store.theta=TRUE, store.alpha=TRUE,
                                     seed = 99813 )

save(out.1d.fixedalpha,
     file="out-1d-fixedalpha-nofraud-balanced-newIDs.Rda")

